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Wc extract scattering phases for fermion-fermion scattering from Monte Carlo simulations of the two- 
dimensional Gross-Neveu model. This is done by means of Liischer's method, which exploits the volume de- 
pendence of the energies of two-particle states. The results are compared with the analytical predictions. 



1. INTRODUCTION 

Calculating scattering phase shifts nonpertur- 
batively is an important problem, especially in 
QCD. Liischer has developed a method to extract 
elastic scattering phases in massive quantum field 
theories from two-particle energy levels in a fi- 
nite volume [1]. Since tlie latter are accessible to 
Monte Carlo calculations, it becomes possible to 
compute scattering phases from numerical simu- 
lations. The procedure has been used successfully 
to study boson-boson scattering in several scalar 
theories (see e.g. [2] and references therein). Here 
we report on the first application to a fermionic 
model: We study fermion-fermion scattering in 
the two-dimensional Gross-Neveu model. 

2. THE GROSS-NEVEU MODEL 

The action of the Gross-Neveu model for n/ 
fiavours [3] reads in the Euclidean continuum 

Sc = j (fx (^alixdn'tpa " yl^aV^a)^^ (1) 

where a = l,2,...,n|. The theory enjoys 
an 0{2nf) symmetry encompassing the obvious 
U(n/) symmetry and charge conjugation. Fur- 
thermore it possesses a discrete chiral invariance, 
which will be called 7,5 in the following. Since this 
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invariance is spontaneously broken, the model de- 
scribes nf massive fermions. The one-particle 
states (fermions and antifermions) transform ac- 
cording to the vector representation of 0{2nf), 
whereas the two-particle states are classified ac- 
cording to the representations of 0(2n/) on ten- 
sors of rank two: the trivial representation on 
invariant tensors (abbreviated as inv. in the fol- 
lowing) , the representations on traceless symmet- 
ric tensors (sym.) and on antisymmetric tensors 
(anti.). The phase shifts in these three symmetry 
sectors are analytically known [4,5]. 

On the lattice we use N copies of stag- 
gered fermions with "hypercubic" interaction [6] . 
Hence the action reads 

- Xa{x + (j.)Xa(x)) 

-f^E (2) 

X \ p ) 

where a = 1,2,...,N and rii,{x) are the usual 
sign factors for staggered fermions. The sum over 
p runs over the corners of an elementary pla- 
quette. For the simulation we rewrite the interac- 
tion in a bilinear form by introducing an auxiliary 
scalar field. Then the hybrid Monte Carlo algo- 
rithm is applied. We work on _L x T lattices with 
the spatial extent L and the temporal extent T 
varying between 8 and 64. 
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The discrete chiral invariance of the lattice ac- 
tion is spontaneously broken, and in the contin- 
uum limit we obtain n/ = 27V massive fermions. 
Again, the continuous flavour symmetry is larger 
than what is immediately apparent. The obvi- 
ous U(A'^) is embedded in the group of unitary 
symplectic transformations in 2N dimensions, 
USp(A'^), which in turn becomes a subgroup of 
the continuum flavour symmetry 0(2n/)=0(4A'^) 
in the classical continuum limit [6]. 

Energy eigenstates with vanishing total mo- 
mentum, which constitute the starting point of 
Liischer's procedure, are classified according to 
irreducible representations of the group of sym- 
metry transformations leaving the time slices 
X2 = const, fixed. The lattice time slice group 
LTS is generated by the USp(7V) transformations, 
the spatial shift, the spatial inversion, and the dis- 
crete chiral transformation. The continuum time 
slice group CTS is generated by the 0(4iV) trans- 
formations, parity, and 75. Note that we can ig- 
nore translations as we are only concerned with 
states of vanishing total momentum. 

The irreducible representations of LTS relevant 
for this work wfll be denoted by A^"""'"'. The 
signs (7i,ai, correspond to the spatial shift, the 
spatial inversion, and the discrete chiral trans- 
formation, respectively; D denotes an irreducible 
representation of USp(A'^) on tensors of rank two 
and takes the values inv. (trivial representation 
on invariant tensors), anti. (representation on 
antisymmetric tensors with vanishing symplectic 
trace), and sym. (representation on symmetric 
tensors). Similarly we have irreducible represen- 
tations A'^'^^ of CTS with signs ap,a^ corre- 
sponding to parity and 75, respectively, and I) 
denoting an irreducible representation of 0(4 A'') 
on tensors of rank two. 

Restriction of an irreducible CTS representa- 
tion to the subgroup LTS will in general lead to 
a reducible representation. It turns out that ev- 
ery A'^j^"^"'^ couples to two continuum symmetry 
sectors K''^"'' with different values of ctp (see Ta- 
ble 1). This is analogous to the appearance of the 
parity partners in QCD with staggered fermions. 
Whenever a given continuum symmetry sector 
can be reached by lattice operators from several 



Table 1: Connection between the irreps of CTS 
and LTS. "Suppressed" contributions are shown 
in italics, "preferred" contributions are under- 
lined. 
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lattice symmetry sectors one has the possibflity 
to check the influence of lattice artifacts. 

3. TWO-PARTICLE ENERGIES 

The needed two-particle energies are extracted 
from correlation functions 

Cij(T) = {Oi{T)O,(0)) , i,j = 1,2,..., r, (3) 

where the operators Oj have the following prop- 
erties: 

• Oi{T) fives on the two adjacent time slices 

X2 = T,T + 1; 

• Oi{T) is bilinear in the fermion fields and 
besides XX terms it may also contain terms 
of the form xX ^-nd xX'i 

• Oi{T) has vanishing total momentum; 
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Figure 1. Generalized eigenvalues for the LTS 
irrep A+~7 on a 16 x 32 lattice. The crosses 
(dots) correspond to positive (negative) values of 

• Oi{t), . . . , Orir) belong to one definite ir- 
reducible representation A^jj"'"' of LTS. 

From the transfer matrix formalism we ob- 
tain the following representation of the correla- 
tion functions Cjj(T) as T — > oo: 

n 

X (e-^"-+e-^"(^--)) . (4) 

Therefore the continuum quantum number trj^^ 
and hence the complete continuum symmetry sec- 
tor to which a particular two-particle energy level 
En belongs can be read off from the oscillating 
or non-oscillating character of the corresponding 
contribution to the correlation function and the 
lattice quantum numbers of the operators used 
((T/, in particular). 

However, in the data (for nj = 8) cither all 
contributions are oscillating or they are all non- 
oscillating or one hardly sees any signal. How can 
this behaviour be explained? First, one observes 
that certain combinations of continuum quantum 
numbers cannot be realized by two-particle states 
due to the Pauli principle: For D = anti. con- 
tributions with ap = +1 are suppressed as is 
(Tp = —1 for D = inv., sym. Secondly, in the con- 
tinuum limit our operators, which are nonlocal in 
time, will couple predominantly to one value of 
(Tp, an effect that was also observed for nonlocal- 



in-time meson operators in QCD with staggered 
fermions. For D = anti., inv. contributions with 
ap = —ai=a[^ '^0/2 preferred, whereas for 
D = sym. the preferred value is ap = a^a[^^''^^^'^ . 
In the end it turns out that each lattice operator 
couples (effectively) to two-particle states from 
at most one continuum symmetry sector (see Ta- 
ble 1). 

We are now ready to discuss the computation 
of the two-particle energies. Since in practice T 
is not extremely large, there are additional con- 
tributions of the form (— l)'^ay + bij in Cij{T). 
They should be of the order e^™f^ and result 
from elementary fermions of mass to/ travelling 
once around the lattice in time direction. We 
eliminate them by subtracting Cjj(r -I- 2). Ex- 
ploiting the symmetry in r we arrive at the r x r 
correlation function matrix 

Cij{r) = Cijir)-Cij{r + 2) 

+ Cij{T-T)-Cij{T-T-2) 

= ^4")A(r,ii;„) (5) 

n 

with 

A{t,E„) = 8e-^"^/2 (-<7,4"))%inh(-ii;„) 
xsinh(E„(T + l-iT)) . (6) 

Next we solve the generalized eigenvalue prob- 
lem [7] 

{v = l,2,...,r) for fixed to {to = 1 in the 
following). From the corresponding eigenvalues 
A^(t) = A(T,E^)/A{to,E^) the energies Ei < 
E2 < ... < Er are then obtained by a one- 
parameter fit. (The influence of the neglected 
states above Ej. can be tested by varying r, the 
number of operators considered.) 

In Fig. 1 we plot In |A^| versus t — to for one 
lattice symmetry sector. Positive values of A^ 
are indicated by plus signs, negative values by 
dots. The curves result from one-parameter fits. 
The number of continuum flavours is 8, and the 
coupling is such that to/ = 0.401(2). 
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Figure 2. Scattering phase from the LTS irreps 
A+JJ+ and AjJ"jj"[j7 (which all couple to the sym. 
sector in the continuum limit, see Table 1) for 
lattices 8 x 8, 16 x 32, 32 x 32. 

4. THE SCATTERING PHASE SHIFTS 

According to Liischer, the scattering phase 6 
at the momentum k^, related to the two-particle 
energy Ei, through 



is given by 



2 " 



modTT 



(8) 



(9) 



(provided polarization effects are negligible) [1,7]. 
So here the continuum dispersion relation enters. 
Since in the fermion propagator the lattice dis- 
persion relation for free formions turns out to be 
very well satisfied we feel encouraged to replace 
the continuum dispersion relation by the lattice 
dispersion relation in an attempt to compensate 
cutoff effects at least partially [2]. Hence we de- 
termine k,, from 



sinh(i£y) = \Jm?f+ sin^ k^ 



(10) 



The scattering phase 5 at is then taken from 
Eq. (9). In Figs. 2,3 we show the results for two 
continuum symmetry sectors coming from a vari- 
ety of lattice symmetry sectors and lattice sizes. 
Thus the plots give an idea of the systematic er- 
rors involved. 

The curves represent the analytically calcu- 
lated phase shifts [4,5]. We observe good agree- 
ment up to k/rrif « 3 showing that Liischer's 
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Figure 3. Scattering phase from the LTS irreps 

^anliT) ^anti^) KntL and A+-7 (wllicll all COU- 

ple to the anti. sector in the continuum limit, see 
Table 1) for lattices 8 x 8, 16 x 32, 32 x 32. 

method works also for fermions - at least in two 
dimensions. The deviations for larger momenta 
may be due the fact that the higher energy lev- 
els do no longer correspond to pure two-particle 
states. 
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